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We use a simplified self-consistent method to address nonlinear-optical cascading phenomena, 
which shows added microscopic cascading contributions in high-ordered nonlinear susceptibilities 
through fifth order. These cascading terms in the microscopic regime encompass all possible scalar 
cascading configurations. The imposition of geometric constraints further influences the predicted 
cascading contributions and opens up additional design parameters for nonlinear-optical materials. 
These results are used in approximating the effective fifth-order susceptibility in thin films of Ceo 
monomers of varying thickness and concentration. This paper contains the corrections to the original 
paper that appeared in the Journal of Physics B as reflected in the content of the corrigendum that 
followed. 



INTRODUCTION 

Materials with large nonlinear-optical responses are 
used for creating devices for harmonic generation, spon- 
taneous down conversion, the optical Kerr effect, etc.[lH3] 
Cascading (at the molecular level) has been a recent topic 
of interest due to the apparent increase in higher-order 
nonlinear optical susceptibilities that stem from products 
of lower-order nonlinearities.[4'-T 

Dolgaleva, et al.^ showed that local field correc- 
tions based on Bloembergen's method [8 can pre- 
dict trends in the nonlinear susceptibilities as functions 
of concentration. [9^ Here, we use a self-consistent 
method to approximate the local field factors and the 
cascading contribution of the first four hyperpolarizabil- 
ities of a geometricly-constrained system of molecules. 
The self-consistent method has previously been shown 
to give the exact solution for the second hyperpolariz- 
ability (third polarizability) in a system of two inter- 
acting dipolar molecules. JTJ [12] There are many advan- 
tages of manipulating the nonlinear-optical response mi- 
croscopically, such as prolonging the stability of optical 
solitons,[13 or adding new design parameters to solid 
film optical limiters.fTT To this end, we introduce a self- 
consistent theory of microscopic cascading for a system 
with a large number of molecules subject to simplifying 
constraints to focus attention on only the most relevant 
tensor components in the dipole approximation. This 
approach is well suited for a gas-lattice model, which we 
use to reduce the translational degrees of freedom when 
calculating large-scale molecular interactions. [ISHTTl 

We then describe the effects of a confined system 
geometry in order to highlight key design parame- 
ters for materials composed of molecules that have a 
large nonlinear-optical response. By increasing the in- 



cident surface-to- volume ratio of thin samples containing 
molecules with a large linear- and nonlinear-optical re- 
sponse, higher-ordered responses can be enhanced via 
an increase in the linear local field and microscopic 
cascading. With recent advances in layered polymeric 
systems, [THHH] multi-layered materials could be created 
with enhanced nonlinear-optical properties due to sur- 
face effects caused by interleaving materials with small 
polarizabilities placed between each active layer. 



THEORY 
Local field effects 

Consider a system of dipoles that are polarized along 
one molecular axis in the direction of an applied electric 
field, which is taken to be the z-axis. The components of 
the molecular susceptibility in the direction perpendicu- 
lar to the polarizable axis are assumed to be negligible. 
In this system, an induced dipole, p^, is polarized by the 
applied field, E^, and the induced electric field of the 
surrounding molecules, X^Ej. 

The problem becomes complicated with five degrees 
of freedom, two rotational and three translational. For 
simplicity, the molecules are assumed to be located at 
points on a cubic lattice. [HI [12] When the only inter- 
action is from nearest neighbors (nn), the electric fields 
inside the lattice from neighboring molecules are simply 
E± = —p/r^ from a neighboring molecule in the direc- 
tion perpendicular to the applied field, and = 2p/r^ 
from a neighboring molecule in the direction parallel to 
the applied field. 
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More interactions can be included in the same fashion 
for second nn, third nn, and so forth. The dipole field in 
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vector form, 
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is used to describe the field from neighboring molecules 
that are not translated either perpendicular or parallel 
to the molecule with respect to the applied field. [22 

The magnitude of the lattice local field factor for two 
interacting molecules, Liat,i, has previously been derived, 
and given by 



^lat,l 
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where /i^2 is a geometric dependent coefficient for two 
identical molecules, a is the polarizability, and r is the 
distance of separation.^! Note that r = |r2 — ri| 
because pi and p2 are nearest neighbors. Equation |2] 
does not include the Lorentz local field, and was derived 
from the self-consistent equation of two identical inter- 
acting dipoles subject to an applied field as opposed to 
the macroscopic field such that 
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FIG. 1: The interaction coefficient, of a molecule at 

the center of the surface of a film as a function of square 
side length for various film thicknesses. The inset shows the 
cubic lattice structure, where the dark striped molecule is the 
molecule of interest. 
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/c^^^ (^Ea + /i,2^^ , (3) which gives the interaction coefficient of the ith molecule. 



where /c^^^ = a and pi = P2 = Pi- 

For an infinite isotropic volume subject to a macro- 
scopic field, Em^ a spherical molecule's direct electric 
field in a cubic lattice is in the direction of Em such that 
the direct field, £^dir, was shown by Lorentz [23] to be 
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where r is the lattice constant and V = Nr^ is the vol- 
ume with N denoting the number of induced dipoles. 
When aberrations in spherical symmetry are present or 
a simple cubic lattice structure is not an accurate physi- 
cal description, then the direct field is defined as 
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where fi^i is the direct field coefficient. Note that when 
defining the polarizability and susceptibility in terms of 
the applied field, we may sum the fields produced by 
all other molecule, which includes the contribution from 
the depolarization field and £^dir in the macroscopic ap- 
proach. 

Assuming that each dipole in a many-body system can 
be approximated as equal and subject to an applied field, 
we can generalize Equation [3] to give 
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Figure [l] shows the calculated f^^ for the gas-lattice 
model in the simple case where the zth molecule is located 
at the center of the surface of a square film. Though lim- 
ited to this single location, the plot illustrates the rate of 
convergence of the interaction coefficient as a function of 
side length for a molecule at the surface and shows that 
this rate depends on the material's ratio of thickness to 
side length. Figure [2] illustrates how of a cen- 

tered non-surface molecule varies as a function of depth 
through a 61-molecule-thick square film for squares with 
various side lengths. 

We define the effective (hyper)polarizabilities, ^ in 
terms of an expansion in the applied field, Ea- These 
(hyper)polarizabilities contain terms that are less than 
or equal to the order in the power series for molecules 
with no permanent dipole moment, [11] as given by Equa- 
tion [3] Therefore, if we assume molecules with no per- 
manent dipole moment, k^^^ ^ 0, and choose to solve 
only for the linear polarizability, agff = then all 

higher-order terms will not be included in the solution. 
We can obtain a solution to the nth-order effective (hy- 
per) polarizabilities by solving for the dipole moment in 
Equation [3j and then make the substitution 



where the effective polarizability for each molecule in the 
cubic lattice with N lattice points is calculated by sum- 
ming over the interactions of all other molecules, j ^ 
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Thus, the effective linear polarizability of the ith 
molecule in a system of N molecules that are subject 
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FIG. 2: The interaction coefficient, f^^ ^\ as a function of 
deptfi in a 61-molecule-thick film for various side lengths of 
square films. Note that the molecules are constrained to po- 
larize only along an axis that is parallel to the applied field. 



to an applied electric field is 

<^eff,i = OiLi, 

where the net local field factor is given by 
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The volume of a primitive cell, Vp = r^^ is the cube of 
the distance between two lattice points in our cubic gas- 
lattice model and summing up all cells gives V = NVp^ 
where V is the volume. Thus, we can rewrite the local 
field factor as 
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where x^^^ is the "undressed" linear susceptibility, x 
Na/V .{x^^'^L is commonly known as the "dressed" linear 
susceptibility when Li = Lj = L for all N molecules. The 
denominator in Equation [TT] can be expanded such that 
Li ^ 1+/^^^ ^^X^^^+' ' ' 7 where the expansion is typically 



truncated to first order since 

yp-i)^(i) ^ ;l ) Note that 
the local field factor in the gas-lattice model is equal to 
the Lorenz field factor when = and written in 

terms of the macroscopic field for a large sphere, where 
the depolarization field is equal and opposite of the direct 
field. This is also observed in the center of a cube of 
61 X 61 X 61 lattice points illustrated by the solid line in 
Figure [2] 



Microscopic cascading 

So far our study has been limited to linear optics. Sev- 
eral groups have reported on the enhancement of third- 
and fifth-order susceptibilities from the microscopic cas- 
cading of lower-order nonlinear interactions. j24]-[28] Al- 
though the molecules in this study are assumed to have 
no permanent dipole moment, an asymmetric potential 
is still possible, which allows even-ordered polarizabilities 
to exist. For example, the hyperpolarizability, (3 = k^'^\ 
can be nonzero without the existence of a permanent 
dipole moment for cases where there exists a nonzero 
octupole moment. [29] 

The nonlinear-optical response of a system is charac- 
terized by taking into account higher orders of n in Equa- 
tion [s] The (hyper) polarizabilities of the iih. molecule, 
which is only allowed to be polarized in one direction 
with all dipole moments approximated as equal in magni- 
tude, can be derived by solving the self-consistent dipole 
moment equation. After algebraic simplifications, the re- 
sulting first five effective (hyper) polarizabilities, ^ are 
given as 



and 



eeff,^ = L'le + 3LjF, (7^ + 2p5) 
+ 21LIf^/3^j^UL^^F^/3^, 

respectively. The local field factor, L^, in Equations [T2 
16 is given by Equation [TT] and Fi is defined as 
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Equations 15 and 16 are given for the first time to our 
knowledge, where all microscopic cascading terms repre- 
sent the possible macroscopic cascading schemes. 

The expressions for the polarizability and first hyper- 
polarizability do not contain any cross-terms. The second 
hyperpolarizability and higher-ordered hyperpolarizabil- 
ities, however, have cross-terms that are ordered by pow- 
ers of Fi. Because Fi is parameterized by the molecular 
interaction coefficient, fj^^ ^\ we find that the cross- 
terms are a fundamental part of a finite volume due to 
the system with boundaries. In other words, when the 
finite volume has a molecular interaction coefficient of 
zero magnitude, the cascading contribution is zero, which 
gives zero contribution from a molecule's hyperpolariz- 
ability to another molecules second hyperpolarizability. 
The most interesting aspect of this calculation is the pos- 
sibility of increasing this dependence by maximizing the 



sum of all interaction coefficients, ^ /• 
ometric constraints. 
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TABLE I: Calculated Molecular Susceptibilities of Ceo 



Symbol 


Value 




Units 


a 


1.85 X 10" 


-23 


cm^ 


7 


3.41 X 10" 


-35 


erg~"'^cm^ 


e 


6.32 X 10" 


-45 


erg~^cm^ 



THE ENHANCED NONLINEAR RESPONSE OF 

Ceo 



The implications of the model in the previous sec- 
tion can be illustrated in a simple analysis of spheri- 
cal molecules. For example, the molecule Ceo is a cen- 
trosymmetric molecule with an approximately spherical 
shape. [30] In the previous sections, a planar geometry 
increased the molecular interaction coefficient, thereby 
strengthening the local field factors and cascading con- 
tributions. Therefore, a thin sheet of Ceo fullerene is 
proposed. Ceo is ideal for our theoretical treatment be- 
cause optical poling is not required to get the maximum 
response due to the molecular configuration, [121 EI] and 
all cross-terms with an even-ordered molecular suscepti- 
bility, such as /3 or (5, are zero from centrosymmetry. 

A thin film of a molecule with a large fifth-order re- 
sponse, such as Ceo, would greatly increase a molecule's 
Li and F^. The nonlinear properties of Ceo have been well 
documented including enhancements to the fifth-order 
nonlinear susceptibility. [32H36] The dependence on the 
third susceptibility squared was shown by Dolgaleva, et 
at. [H [To] Concentration dependence due to aggregation 
was shown to level-off the third susceptibility. [37 Knize 
used Ceo to show how the model of free electrons con- 
fined to a sphere fit experimental results. [38] Moreover, 
it was recently shown that large nonlinear responses were 
obtained from thin films of Ceo -[39] Because "undressed" 
X^^^ effects are much smaller than lower-order material 
responses, one would conclude that the large nonlinear 
response of Ceo in a thin film would be an ideal test case 
to maximize both the local field and microscopic cascad- 
ing contributions. 

The molecular susceptibilities shown in Table [ij are 
found by approximating Ceo as a three level model, [40^ 
[41] where Ceo has two large oscillator strengths with re- 
spect to their corresponding transition energies, which 
are much larger than the oscillator strengths associated 
with other transitions. We stress that the values are es- 
timated from a 3-level model where higher order hyper- 
polarizabilities may be overestimated. The values of the 
transition energies and oscillator strengths were taken 
from the experimental values given by Leach, et a/. [42| [43] 



Interaction coefficient components and fringe effects 

Until now, the self-consistent model derived in the The- 
ory Section has only been used to calculate the effec- 
tive molecular susceptibilities of molecules that are con- 
strained to polarize along a single direction. Ceo, how- 
ever, is a spherical molecule that will equally polarize in 
three dimensions. With this, we define the interaction co- 
efficient in three dimensions, where ^\ ^\ and 

are the respective interaction coefficients in the x, 
and z directions. Although the interaction coefficients 
are expected to be small in the directions perpendicular 
to the electric field, they may not be negligible in certain 
geometric configurations of molecules with large values 
of a/r^. 

Previously, we denoted as the interaction coeffi- 

cient. However, this is the interaction coefficient when all 
induced dipole moments and all interaction coefficients 
are equal in both direction and magnitude. Because the 
interaction coefficient changes as a function of location 
on the lattice, the induced dipole will also change with 
respect to a molecule's position on the lattice. We will de- 
note these interactions as second-order interactions (the 
induced dipole from all other dipole fields that occur 
from first-order interactions). An example is the case 
where the dipole moment of one molecule is subject to 
the dipole fields of the other molecules but not subject 
to the applied field, and the second-order interaction is 
when the dipole fields from the molecules outside the ap- 
plied field contribute to the induced dipole moment of 
the molecules inside the applied field. 

Because the largest field in a lattice of Ceo niolecules 
is produced by the linear polarizability with a magni- 
tude given in Table [T] and the distance of closest ap- 
proach, Tmin, is approximately 7.1 A, [44] we find that 
^/^min ^ 0.05. Therefore, all second-order interactions 
will have a contribution that scales as o? jr^ ~ 2.5 x 10~^. 
Thus we choose to neglect second-order interactions in 
the Ceo lattice based on the small o? jr^ dependence. To 
this end, we focus on the nonlinear susceptibility from 
only the applied field and first-order interactions, and 
assume that ~ through the region illuminated by 
a beam with constant magnitude in that region. Inas- 
much as the applied field is to be taken along the z-axis, 
"Pi^z is the only non-zero component of the induced dipole 
moment produced by the applied field. Due to the possi- 
bility of a anisotropic system geometry, however, all three 
components of may exist due to orientations of a dipole 
field from a molecule in the system with respect to other 
molecules. Although small enough for a valid first-order 
interaction approximation, components of the interaction 
coefficient in the directions perpendicular to the applied 
field are larger near the surfaces of illuminated regions. 
This fringe behavior is illustrated in the blowup of the 
vector diagram of the interaction coefficients shown in 
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FIG. 3: A vector diagram of the three component interaction 
coefficient. The interaction coefficients are the results of an 
off-resonant top hat beam with a cross section of 41 molecules 
along the x- and y-axes, where the beam is illuminating a 7- 
molecule thick film. 



Figure [3] 



Thin films subject to an ofF-resonant top hat beam 



gregation effects. With this assumption, the two exper- 
imentally controllable parameters are the concentration 
and geometry. For this statistical model, we wish to look 
at the fifth-order susceptibility of a thin film of Ceo and 
vary the concentration and thickness. For the sake of 
comparative analysis, the coherent beam's diameter is 
set to be a function of concentration so that there is al- 
ways a constant area of molecules being illuminated. 

The effective fifth-order susceptibility along the direc- 
tion of the electric field for molecules that possess a 
nonzero value for both the second and fourth hyperpo- 
larizabilities is 



N 



= ^E^'+3^^E^'^^ (18) 
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where gi fj;^ 



-1) 



18 



are 



The summations in Equation 
functions of the concentration and geometry. Therefore, 
when ignoring energy perturbations, the effective fifth- 
order susceptibility can be maximized by adjusting these 
two parameters. Although small, it should also be noted 
that the fringe effects discussed in the previous subsection 
become apparent in this geometry, and will give rise to 

(5) 

nonzero components of xli^ ijkimn ^^^^ would otherwise 
be zero. 

Figure W shows the calculated ^ ^^^^ com- 

posed of Ceo molecules subject to a normally-incident 
top hat beam with a frequency far below resonance, i.e., 
dEa/dx ~ 0. Here, xlff is calculated as a function of 
concentration for different film thicknesses. 

The effective fifth-order susceptibility shown in Figure 
|4]is calculated by summing the effective fourth hyperpo- 
larizabilies of every molecule in the top hat beam radius 
and dividing by the cylindrical volume. Although the 
number of molecules in the calculation is small relative 
to large collimated beams, one can predict the values of 
Xeff for thin samples in the macroscopic regime by scal- 
ing the fraction of thickness to beam diameter, thereby 
strengthening the approximation of a negligible fringe ef- 
fect near the beam edge. 



We stress that most of the terms in Equation [16] do 
not contribute to the fifth-order susceptibility of Ceo due 
to centrosymmetry. One would be correct in concluding 
that an asymmetric molecule would greatly increase the 
number cascading terms, although this does not neces- 
sarily mean that the total fifth-order nonlinear response 
will always be larger. 



The effective fourth hyperpolarizability, Ceff, of the cen- 
trosymmetric molecule Ceo depends on the polarizabil- 
ity, second hyperpolarizability, fourth hyperpolarizabil- 
ity, Ceo concentration, and geometric constraints. The 
polarizability, second hyperpolarizability, and fourth hy- 
perpolarizability are assumed to be constants of Ceo- 
This is true so long as we ignore perturbations of the 
energy levels due to molecular interactions such as ag- 



Approximations and Assumptions 

The results of these calculations are based on a num- 
ber of approximations. The model does not account for 
aggregation effects since each molecule is required to be 
positioned on a lattice point at a fixed separation. There- 
fore, large energy perturbations and the potential for 



6 



many molecular systems to "level-off" the nonlinear re- 
sponse at high concentrations are not considered. A com- 
plete theory of cascading contributions and the influence 
of molecular interactions on the transition energies has 
been previously derived for the two molecule system, [12] 
which can be folded into the current model with some 
additional rigor. A full tensorial treatment up to third 
order of the cascaded nonlinear effects with molecular 
correlations is in ref [4 . 

An off-resonant response is assumed such that the fre- 
quency of the electric fleld approaches zero. Thus, retar- 
dation effects over large distances as well as transitional 
behaviors near resonance are not treated. Also, because 
there are fleld gradients, albeit small, there will be ef- 
fects from the dipole-quadrupole polarizability, dipole- 
octupole polarizability, quadrupole-quadrupole polariz- 
ability, etc., leading to deviations from the dipole ap- 
proximation. Finally, all dipole moments are taken as 
equal to allow a self-consistent approach. This assump- 
tion ignores higher order interactions, which may play an 
important role for asymmetric systems that contain very 
small numbers of molecules. 

As an additional point, the electric fleld is assumed to 
be uniform throughout a cylindrical volume, but vary- 
ing the magnitude of the applied fleld over the substrate, 
such as a thin fllm subject to various beam proflles, would 
allow one to approximate the cascading contribution of 
an optical material for a multitude of experimental con- 
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ditions. Under these circumstances, the dipole moments 
can no longer be assumed to be equal and other assump- 
tions must be made about the spatial distribution of 
dipole moments induced by the applied fleld. For ex- 
ample, a Gaussian beam proflle can be approximated by 
a Gaussian distribution of dipole moment magnitudes. 



CONCLUSION 

The self-consistent method used to evaluate the ef- 
fective nonlinear susceptibilities provides a solution that 
reveals multiplicative factors that are the same as the 
Bloembergen's approach provided that the power series 
representation of nonlinear optics is consistent. These 
multiplicative factors show that the cascading contribu- 
tion is not just the sum of lower order terms multiplied by 
a local fleld factor, but that some cascading contributions 
are larger than others. Furthermore, the microscopic cas- 
cading terms that arise from the self-consistent method 
encompass all possible macroscopic cascading schemes, 
and are ordered in terms of F^. 

We have shown that it is possible to increase the lo- 
cal fleld factors and cascading contributions in higher- 
ordered nonlinear responses by forced geometric con- 
straints. Inasmuch as the local fleld factor can increase 
the bulk nonlinear response, a decrease in the distance 
of a molecule from a surface that is perpendicular to 
the direction of the applied fleld can greatly increase a 
molecule's interaction coefficient with respect to the sur- 
rounding molecules. 

The molecule Ceo fullerene was used as a case study, 
which showed an increase in the effective flfth-order sus- 
ceptibility via molecular interactions. Moreover, this en- 
hancement from molecular interactions was shown to in- 
crease with a decrease in molecular separation distances. 
Thus, increasing the concentration will not only increase 
the nonlinear response of the polarization density based 
on number density, but also will show an increase in 
nonlinearity due to interactions at relatively short dis- 
tances. These interactions are shown to further increase 
due to geometric constraints imposed on the active vol- 
ume. The model explicitly shows how a material can 
achieve a larger effective nonlinear response than pre- 
dicted for a bulk material. 
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